Isolated large amplitude periodic motions of towed rigid wheels 
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Abstract 



This study investigates a low degree-of-freedom (DoF) mechanical model of shimmying 

wheels. The model is studied using bifurcation theory and numerical continuation. Self-excited 

vibrations, that is, stable and unstable periodic motions of the wheel, are detected with the 

help of Hopf bifurcation calculations. These oscillations are then followed over a large parame- 

^*^ . ter range for different damping values by means of the software package AUT097. For certain 

r \ ' parameter regions, the branches representing large amplitude stable and unstable periodic mo- 

'^ , tions become isolated following an isola birth. These regions are extremely dangerous from an 

P^ ■ engineering view-point if they are not identified and avoided at the design stage. 

1 Introduction 

►^ , Shimmy is a common name for the lateral vibration of a towed wheel. This has been a well-known 

QQ ' phenomenon in vehicle systems dynamics for several decades: the name shimmy comes from a dance 

CN , that was popular in the 1930s. One of the early scientific studies of shimmy dates back to this time 

^N ' [I]- In many cases, the appearance of shimmy is a serious problem, for example, in case of nose gears 

of airplanes or front wheels of motorcycles. There are many mechanical models (see, for example, 
[H [SI in [11 [51 [7] ) that describe the shimmy of rolling wheels. The two most important considerations 
of any model involve whether the wheel is rigid or elastic and whether the suspension system is rigid 
or elastic. The simplest combination of a rigid wheel and a rigid suspension system does not, by 
definition, support lateral vibrations. Lateral vibration can occur in models with an elastic wheel 
and a rigid suspension (see [8j) and in models where both the wheel and the suspension system are 
elastic. But in this paper we wish to make analytic progress in order to obtain a clear understanding 
$H ' of the dynamics and so we consider a low degree of freedom (DoF) mechanical model of a rigid 

wheel, which has a viscously damped elastic suspension. This model was studied without damping 
in [3], where subcritical Hopf bifurcations and chaotic and transient chaotic oscillations were found. 
This paper is structured as follows. First, in Section 2, the mechanical model is introduced, 
together with the equations of motion. In Sections 3 and 4, the Hopf bifurcation calculation is 
presented in the presence of viscous damping at the suspension. The effect of damping on the 
stability of stationary rolling is analysed. In Section 5, the periodic solutions of the system are 
followed using AUT097 lOj, the stability charts and the bifurcation diagrams are plotted and also 
compared to available analytical results. 

The bifurcation diagrams show an isola birth where isolated large amplitude stable and unstable 
periodic motions coexist with the stable stationary rolling solution. These motions are difficult to 
detect either by numerical simulation or by conventional stability and bifurcation analysis. The 
presence of unstable periodic motions indicates a dangerous system configuration. 

2 Mechanical model 

The mechanical model under consideration is shown in Figure [1] The plane of the rigid wheel is 
always vertical to the ground, and they have a single contact point at P. The radius of the wheel is 
R, the mass of the wheel is m^ and its mass moment of inertia with respect to the z axis is J^z ■ 



The mass moment of inertia with respect to the y axis of its rotation is Jwy^ where subscripts w 
refers to wheeL The caster length is I, the distance of the center of gravity C of the caster and the 
king pin at A is Ic- The mass of the caster is rric and the mass moment of inertia with respect to 
the z axis at C is Jcz, with subscripts c referring to caster. The system is towed in the horizontal 
plane with constant velocity v. The king pin is supported by lateral springs of overall stiffness k 
and viscous damping coefficient b. 

Without rolling contraints, the system has 3 degrees of freedom, so one can choose the caster 
angle '0, the king pin lateral position g, and the wheel rotation angle Lp as general coordinates. The 
constraint of rolling (without sliding) means that the contact point P has zero velocity. This rolling 
condition leads to two scalar kinematical constraint equations in the form of coupled first order 
nonlinear ordinary differential equations (ODEs) with respect to the general coordinates. 

The equations of motion of this rheonomic and non-holonomic system can be derived with the 
help of the Routh-Voss equations or the Appell-Gibbs equations [TT]. In the case of zero damping 
(& = 0), they are given in [12j . Here, we present the equations for non-zero viscous damping: 
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are odd and even functions respectively of the general coordinates. The first two equations in (1) 
are the equations of angular momentum and the last two express the constraint of rolling (without 
sliding) of P. 

Since the general coordinate ip appears only in the fourth equation of motion, this coordinate 
is a so-called cyclic one and the system can be described uniquely in the three dimensional phase 
space of the caster angle ip, caster angular velocity i? and the king pin lateral displacement q. 

3 Stability analysis 

The trivial solution of the system is stationary rolling along a straight line defined by the vector of 
the towing velocity: 

ip = 0,ip = 0,q = and ip = ^ . 
When the towing speed is zero, the system forms a 1 DoF oscillator about the z axis at P. The 
corresponding angular natural frequency uj^ of the undamped linear system and the damping ratio 
C are given by 



y Jiv2 + Jc + mc{l - IcY ' 2 k 

Let us introduce the new dimensionless parameters: 



kP ^ lb .^. 



i = ^ y = z;7r.^ (5) 



where L and V are the dimensionless caster length and towing speed respectively, and k, and x 
are dimensionless mass moment of inertia parameters related to the caster and wheel geometry and 
inertia. 

With these new parameters, the third order Taylor series expansion of the first three governing 
equations Eq. (1) about the trivial solution assumes the form: 
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The characteristic equation is obtained from the linear coefficient matrix of Eq. (7): 
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According to the Routh-Hurwitz criterion, the stability of the stationary rolling is equivalent to: 



L > LcriV) 



V{1 - 4C^ - 2CKy) 



(9) 



2C + kV 

considering positive parameter values only. At fixed damping ratios C, the stability boundary curves 
are characterized by 



K.xt, = 
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(10) 



where dLcr{V)/dV = at V^ = \4xt and Lcr{V) = at F = V^nax- Stationary rolfing is always 
stable for V > T4iax or L > icr(Kxt)- Figure [Ha shows the corresponding stability chart in the 
plane of the dimensionless towing speed V and caster length L for k = 0.203 and x = 5.67. The 
stability region is shaded for ( = 0.1. The dot-dash line is the locus of 14xt as a funcion of (. Our 
parameters come from a realistic towed wheel of a shopping trolley. The parameters of the wheel are 
m^ = 0.3519 [kg], J^y = 4.63 • 10"^ [kgm^], J^^ = 2.38 • 10"^ [kgm^] and R = 0.04 [m]. The data of 
the caster are nic = 0.0668 [kg], Ic — 0.012 [m] and Jcz = 3.48 • 10^^ [kgm^]. We shall discuss Figure 
|2lb in Section 4. 

The stability chart in dimensional terms is of course of great practical use to engineers. However 
the natural choice of scaling in Eq. (5) makes it difficult to immediately deduce its form from Figure 
[2ja. Hence we redraw this Figure in dimensional terms in Figure [3j using the shopping trolley 
parameters. 



4 Hopf bifurcation 

Despite giving a more useful stability chart, the dimensional form is less amenable to analysis. In 
dimensionless form the stability calculations are simpler and clearer. 

The eigenvalues of the linear coefficient matrix of the dimensionless system can be determined 
on the stability boundary when L = Lcr as given in Eq. (9). There are two complex conjugate 
eigenvalues with zero real part and one negative real eigenvalue: 



Ai,2 = ±i^ , A3 = -cj„ {2( + K,V) , 



where 



Vl - 4(2 - 2CkV ' 



(11) 
(12) 



hence, there is Hopf bifurcation on the stability boundary. The eigenvectors {si, S2, S3} can also be 
determined: 
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The transformation matrix T can be constructed from the eigenvectors in the foUowing way: 

T = [ Resi Imsi S3 ] . (14) 

Let us introduce new variables [xi, X2, x^] such that: 



(15) 



If we substitute this into Eq. (7) and we muhiply the equation on the left with the inverse of 
the transformation matrix, the Poincare normal form is calculated, which has the form 
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(16) 



Since the nonlinearities are symmetric (i.e. there are no second degree terms in the nonlinear 
part of the Poincare normal form), the centre manifold is approximated by a second degree surface. 
Thus, the transformation of the nonlinear part needs only the terms in xi and X2- The terms in 
which X3 appear can be neglected. The sense of the Hopf bifurcation comes from the reduced form 
of Poincare- Lyapunov parameter S for the symmetric case (see 13 J: 
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If 6 is positive/negative then the periodic orbit of the Hopf bifurcation is subcritical/supercritical, 
namely unstable/stable. If we set ^ = 8 and take k ~ 0.203 and x = 5.67 as before, S can be plotted 
as a function of the damping C,, see Figure SI In our case the sign of S changes at a critical value of 
the damping. 

To compare our results with numerical continuation, the theoretical bifurcation branch has to be 
calculated. Accordingly the real part of the implicit derivative of the characteristic equation with 
respect to the bifurcation parameter L has to be determined at the critical parameter value L^r'- 
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The amplitude of the periodic orbit is given by 



ReA'l^ . 
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namely 
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With the help of the transformation matrix T, the amplitude of the vibration can be plotted with 
respect to the general coordinates i/j, ■& and q. The theoretical branches of the amplitude of ip 
(A^) are compared with the numerically continued results of the dimensionless system Eq. (7) in 
Figure [2lb. In this figure, as in the remainder of the paper, dashed lines mean unstable branches 
and continuous lines mean stable solutions, and thick and thin lines show numerical and theoretical 
results, respectively. When the bifurcations are supercritical, the numerical continuation shows the 
existence of folds in certain parameter regions, leading to large amplitude unstable periodic motions. 
As mentioned at the start of this section the analysis of the damped dimensional system is 
difficult. But in case of zero damping the sense of the Hopf bifurcation, the Poincare-Lyapunov 
parameter S, can be determined: 



S^- ;; , ^ , ^ . — ^77^2 ( "^™^ + ^'^c^c + ^^J«,j, ) >0, (22) 
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which is always positive for all real parameter values, namely the Hopf bifurcation is always sub- 
critical (i.e. the periodic orbit is always unstable). The comparison of the theoretical and numerical 
results is shown in Figure [5l 

5 Numerical Continuation 

Returning now to the damped case (with 6 = 0.2 [Ns/m]), the branch of unstable periodic solutions 
was followed using AUT097 [TD] from the Hopf point and a saddle-node bifurcation (fold) was 
detected, see Figure IHl This means that the damped system also has stable periodic solutions. So 
if the system is perturbed enough then it will be attracted to these large amplitude stable periodic 
solutions. 

If the damping factor is increased, the peaks of the stable and unstable branches move closer 
together. At a critical damping factor (6 = 0.269 [Ns/m]) the branches intersect each other and an 
isola is born, see Figure [H Simulated motions are shown in phase space in Figure [T] So a separated 
periodic solution branch occurs in the bifurcation diagram for b > 0.269 [Ns/m]. There are now 
three folds in the figure and, by changing the damping factor b or the towing velocity w, the location 
of these folds also changes. 

In practice such large amplitude stable periodic motions might be expected to slide, once a 
certain critical friction force is reached. Although we consider rolling, slipping is not part of our 
model and this would have to be included if we wished to be certain that such motions could be 
observed in practice. 

For fixed damping ratio a nonlinear stability chart can be plotted in the (v, I) plane, where the 
location of the folds can also be marked. With different projections, the nonlinear behaviour of the 
system is represented in Figure [H 

The bistable area is of great importance. It has been found only with numerical continuation. 
That means that the very dangerous large amplitude periodic motion can not be discovered analyt- 
ically in this case. This bistable area can be reduced and bounded for large enough damping ratio. 
This is shown in Figure [9l 



6 Conclusions 

The importance in the investigation of the assumption of a rigid wheel model can be questioned, 
because a real system is more complex. However the point of our study is to show that the well- 
known properties of shimmy, for example the existence of unstable large amplitude periodic motions, 
are present even in this simple model. 

In this paper the damped model was considered using analytical and numerical methods. It 
was shown that the subcritical Hopf bifurcation can change to supercritical if the damping ratio 
is increased. In the damped system a separated branch of periodic large amplitude solutions was 
detected using numerical continuation. For some parameter values, the stationary rolling solution 
is stable for any value of the towing speed, but there is also a separated periodic solution branch, a 
so-called isola of large amplitude, which is very dangerous. This region can not be determined by 
linear analysis. Therefore, without the nonlinear numerical investigation a system may be designed 
with parameter values that allows for the presence of these dangerous solutions. 
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Figure 1: Model of rigid tyre. 
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Figure 2: The stability chart and the bifurcation branches of the dimcnsionless system 
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Figure 3: Stability boundaries of the linear system in dimensional form using the shopping trolley 
parameters. 
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Figure 4: The S parameter in Eq. (18) as a function of the damping ^. 
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Figure 5: The theoretical and numerical bifurcation branches of the undamped system. 
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Figure 6: Bifurcation diagrams of the damped system with the isola birth. 
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Figure 7: The trajectories in the phase space for b = 0.269 [Ns/m]. 
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Figure 8: StabiHty chart of the nonhnear system. 
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Figure 9: Stability boundaries at different damping factors. 
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